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f length = 64; //fast fourier transform length 



400 

/ 



delta = .1 ; %input ('Enter stop band tolerance: 
beta = 0.2; %input ('Enter pass band tolerance: 
wp = 8; 
ws = 7; 

mu = 28; %CP length 

load channeMmpulse_resp.dat; %load channel impulse response 
N = 19; % sirf filter length 
41 0 { ch = channel_impulse_resp; 
mm = max(abs(ch)); 
ch = 2047*ch./mm; 
NN =-size(channelJmpulse_resp); 
M = NN(1); 
for i = 1:mu 

mv(i) = 1 ; %initialize mv 
end 

C = conv(mv,abs(ch)); 
\^ [Y,l] = max(abs(C)); %find max which is beginning of Gl (CP) 



/ % construct the channel impulse response matrix 
/ fori = 1:M 



forj = 1:N 

if i - j + 1 <= 0 

break; 
else 

H(iJ) = ch(i-j + 1); 
end 
end 
end 
n = 1; 

forl = M + 1:M + N-1 
430 < n = n + 1; 

k = 0; 
for J = n:N 

k = k + 1; 

H(i,j) = ch(M-k + 1); 

end 
end 

for 1 = 1:mu 

mv(i) = 1; 
end 

C = conv(mv,abs(ch)); 
\ [Y,l] = max(abs(C)); 
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for i = 1:mu 
mv(i) = 1 ; 
440 { end 

C = conv(mv,abs(ch)); 
[Y,l] = max(abs(C)); 



V 



% initialize tine SIRF filter to an arbitrary value 
c = [ -0.0593 
0.01047 
-0.062386 
0.02418 
-0.065114 
0.030031 
-0.039083 
0.13789 
-0.10266 
-0.014681 
-0.11224 
0.0546 
-0.12642 
0.18608 
-0.020895 
0.38407 
-0.30117 
-0.37885 
0.42326]; 
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g = fft(c,length); 
f = abs(g); 
450 < lamda = 1; 
tol = 50; 

cold = sum(abs(c)); 
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%begin iteration 

for ii = 1 :20 % 20 is tlie number of iterations 
% projection on the set C2 
for m = wp:length/2 
if ( f(m) > (1 + beta)) 
gg = g(nn); 

g(m) = (1 + beta) /f(m)*complex(real(g(m)),imag(g(m))); 
g(m) = gg + lannda*(g(m) - gg); 
end 

if ( f(m) < (1 - beta)) 
gg = g(m); 

g(m) = (1 - beta)/f(m)*complex(real(g(m)),jnnag(g(m))); 
g(ni) = gg+lamda*(g(m) - gg); 
end 

end 

for m = 1 :ws 

if (f(m) > delta) 
gg = g(m); 

g(m) = delta / f(m)*complex(real(g(m)),innag(g(m))); 
g(nn) = gg+lamda*(g(m) - gg); 
end 

end 

for m = 2:length/2 

g(m + length/2) = conj(g(length/2 - m + 2)); 
end 

or = real(ifft((g),length)); %transform to time domain 
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% time domain projection on C4 set 
fori = 1:N 

c(i) = cr{i); %project into the set C4 

end 

% time domain projection on C5 set 
for n = I - mu - 10:1 - mu 
norm = 0; 
prod = 0; 
%prod1 = 0; 
fori = 1:N 

norm = norm + H(n,i) ^ 2; 

prod = prod + c(i)*H(n,i); 

%prod1 = prod1 + c(i)*ch(n - i + 1 ); ^ 
end 

if (prod > tol) 
for nn = 1:N 
cc = c(nn); 

c(nn) = c(nn) + {(tol - prod /norm) * H(n,nn) 
c(nn) = cc + lamda*(c(nn) - cc); 
end 
end 

if (prod < -tol) 
for nn = 1 :N 
cc = c(nn); 

c(nn) = c(nn) + ((-tol - prod /norm) * H(n,nn 
c(nn) = cc + lamda*(c(nn) - cc); 
end 
end 
end 

for n = l;l + 100 
norm = 0; 
prod = 0; 
for i = 1:N 

norm = norm + H(nJ) ^ 2; 

prod = prod + c(i)*H(n,i); 
end 

if (prod > tol) 
for nn = 1 :N 
cc = c(nn); 

c(nn) = c(nn) + ((tol - prod) /norm) * H(n,nn); 
c(nn) = cc + lamda*(c(nn) - cc); 

end 
end 

if (prod < -tol) 
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for nn = 1:N 
cc = c(nn); 

c(nn) = c(nn) + ((-tol - prod) /norm * H(n,nn); 
c(nn) = cc + lamda*(c(nn) - cc); 
end 
end 

end 

ss = abs(cold - sum(abs(c))) 
cold = sum(abs(c)); 

g = fft(c, length); %transform to frequency domain 
f = abs(g); 



\^ end 



for 1 = 1:N 

480 < sirf(i) = c(i); %sirf will hold the SIRF filter coefficients 
end 



FIG. 4E 




FIG. 5 



